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SUMMARY 

The generalized potential that gives the scattering into a selected group of inelastic 
channels along with the elastic scattering is derived by means of the projection operator 
for the set of channels. A general class of projection operators that selects all open 
channels is developed in terms of the previous projection operator. The resonance part 
of the generalized potential can then be obtained by means of Feshbach's unified theory of 
nuclear reactions. There is no limitation on the types of reactions that can occur. Ex- 
amples of projectors belonging to the class are given for the case of two-channel reac- 
tions. An example which is discussed in some detail is the pickup process. After the 
resonance contributions have been isolated, the transition amplitudes can be energy aver- 
aged. The result is used to obtain the generalized optical potential. The properties of 
the potentials that follow from the equations which determine them are discussed. A 
short review is given of some elements of Feshbach's reaction theory that are required. 

INTRODUCTION 

The complex potential model has proved to be very useful for representing the elas- 
tic scattering of particles by complicated target systems in a wide variety of situations. 
More recently, with the continuing improvement in computer technology, the use of gen- 
eralized potentials has become practical and has given equally impressive results. Many 
derivations of the exact single- channel potentials have been given (refs. 1 to 7). They 
lead to formal expressions for a potential which gives precisely the same elastic scatter- 
ing as that determined by the many-body Schrodinger equation from which the probelm 
starts. 

Multichannel or generalized potentials have been derived by Feshbach (refs. 1 and 2) 
and Mittleman and Pu (refs. 8 and 9). None of these treatments, however, combine all 

*To be published in the Physial Review. 
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the following conditions: 

(1) There is no limitation on the types of reactions that can occur. 

(2) Resonance contributions to the potential are isolated and shown explicitly. 

(3) The potential obtained describes the elastic and any selected part of the inelastic 
scattering. 

In medium- energy nuclear reactions, however, all these conditions must be met in order 
to relate directly an exact potential to a phenomenological potential of the kind employed 
in a coupled- channel calculation for the elastic and inelastic scattering. The primary 
purpose of this report is to show how the generalized potentials can be derived when all 
three of the preceding conditions are fulfilled. In the course of treating this problem, 
several useful properties of channel projection operators are developed. Furthermore, 
a rigorous procedure for including resonances in coupled- channel calculations is exhib- 
ited as part of the overall development. 

The term generalized potential refers to a simplified Hamiltonian which will yield 
a scattering solution with the same asymptotic behavior in a certain group of channels as 
that which the scattering state of the full Hamiltonian, with the corresponding boundary 
conditions, has in these channels. The term generalized optical potential is reserved to 
designate a complex potential that produces, at an energy E, transition amplitudes that 
are equal to the energy-averaged transition amplitudes of the system when the average 
is taken over a suitable energy interval about the value E. (Symbols are defined in 
appendix A. ) The counterpart employed in applications is a phenomenological potential 
model that possesses internal degrees of freedom corresponding to a particular kind of 
target excitation. The model is used to describe the inelastic scattering for a selected 
set of target excitations of the type represented by the particular model, as well as the 
elastic scattering. The approach leads to a set of coupled equations for the relative mo- 
tion part of the wave-function components in these channels. Solving the set of equations 
yields, concurrently, the scattering given by the model in all the included channels. 

In direct nuclear reactions, there frequently occur inelastic scattering amplitudes 
for the excitation of highly collective states of the target nucleus which are almost as 
large as the elastic amplitude. In such cases, approximations for the inelastic scattering 
which neglect the coupling between these strongly excited channels are often inaccurate. 

At the same time, the single- channel optical model may require unusual values for its 
parameters in order to describe the direct elastic scattering. A set of coupled equations 
must be used to obtain satisfactory results. As the findings of coupled- channel calcula- 
tions (refs. 10 to 12) accumulate, it is expected that the empirical properties of general- 
ized optical potentials will become as well established as those of the single- channel op- 
tical model and provide information about the nucleus of comparable importance. 

In the next section, GENERALIZED COMPLEX POTENTIALS, a generalized potential 
is derived (in a way that will be useful later) that gives the scattering in a chosen set of 
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inelastic channels along with the elastic scattering. This result has been given recently 
by Mittleman and Pu (refs. 8 and 9) and it satisfies conditions (1) and (3). This kind of 
potential, however, is not well suited to the situation where the scattering exhibits reso- 
nance effects because of the formation of quasi- stationary bound states of the total sys- 
tem. In this case, what is usually needed for comparison with experimental results is the 
optical potential. To obtain the optical potential, it is first necessary to separate explic- 
itly the resonance part of the transition amplitude, which contains the rapid energy de- 
pendence, from the slowly varying part. For this purpose the formalism of Feshbach's 
unified theory of nuclear reactions (refs. 1 and 2) is particularly convenient. Feshbach 
has obtained expressions for generalized potentials, both the monoenergetic and optical 
variety, which describe the reactions into all the open channels. In reference 1, the de- 
velopment is limited to the case where only inelastic reactions occur and thereby satis- 
fies conditions (2) and (3). In reference 2, the treatment is extended to general reactions 
and so meets the requirements of conditions (1) and (2). In appendix B, some basic ele- 
ments of Feshbach's reaction theory that are required herein are reviewed. 

In this report the potentials which satisfy all three of the preceding conditions are 
derived. The key step in the development is to relate the projection operator P, which 
selects all the open channels, and the projection operators employed in the section 
GENERALIZED COMPLEX POTENTIALS, which select the inelastic channels. This is 
accomplished in the section SUBSPACES AND THEIR PROJECTORS. First, the proper- 
ties of the P projectors are described; a form for P is chosen that has the desired 
relation; and then the general expressions are applied to reactions that involve only two 
types of channels. By way of illustration, the results are discussed for the case of pickup 
or stripping processes in appendix C. The generalized potential, which describes the 
scattering in the selected inelastic channels, is obtained in a form that clearly exhibits 
the resonance contribution. Appropriate energy averaging then leads to the correspond- 
ing generalized optical potential. Finally, some of the properties of these potentials that 
follow from the equations which determine them are discussed. 


GENERALIZED COMPLEX POTENTIALS 

The collision of particles that are themselves bound aggregates of elementary parti- 
cles is now considered. A reaction channel is identified by the indices a and n which 
designate, respectively, the distribution of fundamental particles into bound fragments 
and the internal states of these fragments. For simplicity, all particles are treated as 
distinguishable. Then, for a given set of bound fragments a, the Hamiltonian can be 
separated into a noninteracting Hamiltonian and an interaction V^; 
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( 1 ) 


Jf=K + V 


a 


a 


The operator K a contains the kinetic energy of each elementary particle and its interac- 
tion with every other particle of the bound fragment in which it is found in the set a and 
may conveniently be expressed as 


K a = T or 


+ Jt 
a 


( 2 ) 


where each term on the right consists of a sum over all the a fragments, the sum 
of the center-of-mass kinetic energies, and j? the sum of the internal Hamiltonians. 
An eigenstate (P an of consists of the product of the state of motion of each frag- 
ment relative to its fixed mass center. An eigenstate u ? of can be taken to be a 
product of plane waves corresponding to the free motion of the mass centers, where k is 
used to specify the momenta of the fragments. Noninteracting states, which are eigen- 
states of K^, can be written as 


%nk = u ak%n 


( 3 ) 


The indices ao will be reserved for the incident channel. Of course, incident chan- 
nels consisting of two particles are of primary interest, but the formalism appears the 
same in the more general case. The continuum eigenstate of jf, that describes the re- 
action, is given by the formal expression 


^aok %ok + e - + j£ Va ^ aok ^ 

where the limit e — 0 + is understood on the right. The plus sign on the term ie cor- 
responds to the boundary condition that only outgoing waves are present at asymptotic 
separations of the fragments in the exit channels /3m * ao. The center-of-mass refer- 
ence system with the center of mass at the origin is the most useful reference system 
for our purposes. For notational convenience, the abbreviated form 


* = * 


<+)_ 

a ok 


is used for the desired solution (eq. (4)) of the Schrodinger equation 

(E - jf )* = 0 


( 5 ) 
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In addition to the explicit form of equation (4), the Lippman- Schwinger integral equation 
is satisfied by 


1 * r ^aok 


E-K a+ ie 


V 'k 
a 


( 6 ) 


Consider now the subspace ^ arn of the total Hilbert space spanned by all states of 

the form X<P am , where x is any state of motion for the mass centers of the particles a. 

Because any two <P am with different values of m are orthogonal on the scalar product 

over the internal variables of the particles, the constitute a set of subspaces which 

which are orthogonal on the full scalar product of the space. The projection operator 

n into the subspace tT ^ can be written as the outer product 
am am 

\m' | cum 

where it is understood that n acts as a multiple of the unit operator on a state of mo- 
tion of the mass centers. Then the set of projectors has the following properties: 




(7a) 


n 


t 

am 


= n 


am 


(7b) 


^am^am 1 ? ^mm^am'f 


(7c) 


and 


[ n am' K J = ° <7d) 

The wave function can be expanded in terms of the complete orthonormal set of eigen- 
states of 


\jr = \ f (p (8) 

w / v am^am v ’ 

m 

where the coefficients f are functions of the center- of- mass positions of the particles 
in the set a. The sum runs over the complete set of eigenstates <P am including the 
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states in the continuum spectra of It then follows directly that 

" Warn ^ 

and that all the information about the asymptotic behavior of & in the channel am is 
contained in the projection II Q!m 1 J r . 

A set I of a channels, which includes the incident channel ao but otherwise con- 
sists of an arbitrary selection of a channels, is now considered. The operator 

n l - E n am < 10 > 

mcl 

projects into the vector or direct sum of the selected orthogonal spaces: 

r i = r ao® r ai® • ' ' ®*aj °» * • •» 3 cI 
The projector into the entire subspace orthogonal to ^ is 

©j= i - n r Ojnj= o (ii) 

The projections of ^ are 

<p T - iir 1 ^ = S f cp 
I I / y am r am 

mcl 


and 


= 0j^ 

Operating on equation (5) with IIj and 0j gives the pair of coupled equations: 

(E - (12a) 


and 


(E - Oj-JfOj)^' = 


(12b) 
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From the commutation relation (eq. (7d)) it is seen that 


n^iij = K^iij + njV a n r n^Qj = njV a o p etc. 


Equation (12b) can be solved as follows: 


** = 0 T V„* 

E - K a + fe L 01 


( 13 ) 


This solution is just the result obtained by operating on equation (6) with 0j. The plus 
sign on the ie term ensures that, at large separations, only outgoing waves appear in 
the open inelastic channels not contained in I and in the contribution of 4'* in the re- 
arrangement channels. The imagninary term can be omitted if I contains all the open 
channels. The vector V Q ,\f r satisfies the equation 


V ^ = V + V 
a a I a 


E 


K a+ ie 


0 I V «* 


If an operator 



is defined by the equation 



= Vo. 11 ! + V « 


E - K a + ie 


0 I U I +) 


n. 


(14) 


(15) 


it is clear that 


= uJ + ^j 

satisfies the same equation (eq. (14)) as does V Q ,'i r and, hence, that 

V a* = U I +) *I ( 16 ) 


The propagator in equation (14), 

(E - K a + ie) -1 ©! = 0j(E - K a + fe) -1 ©! 


can be rewritten as 
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(E - K a + ie)' 1 ©! = ©j(E - K a - ©^ + fe) -1 ©^ - V a (E - + fe)' 1 ©^ 

= (E - 0j^+ ie)" 1 ©^ - V a (E - K a + ie)' 1 ©^ 

Inserting the resultant identity into equation (15) gives an explicit formal solution of the 
equation, namely. 


U 


!+) = ww 


©J^ + ie 


0 I V <* n i 


(17a) 


This last result could have been obtained directly by noting that the solution of equa- 
tion (12b) can also be written as 

= (E - ©jjfQj + k)" 1 ©^^ 

The generalized potential is defined as 

V I +) = n i U l +) (17b) 

and has the property that 

"iV - 

Equation (12a) can then be written in the form 

( E - K a- vj +, )*i=° (18) 

which shows that the potential will give the true scattering in the channels I. The 

generalized potential of equation (17) has recently been applied by Pu (ref. 13) to the 
scattering of electrons by atomic hydrogen. 
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SUBSPACES AND THEIR PROJECTORS 
Properties of Projectors That Select All Open Channels 

The generalized potential given by equation (17) is not convenient when the scattering 
shows resonance behavior. The potential will then have a rapid energy variation corre- 
sponding to poles in its energy dependence. The origin of this effect lies in the possibility 
of the system forming compound states with an appreciable lifetime for decay into the exit 
channels. In order to isolate the contribution to a generalized potential arising from this 
aspect of system behavior, one of the fundamental theories of reactions must be employed, 
which will provide a formalism for treating resonances. The theory of nuclear reactions 
due to Feshbach is particularly appropriate, because it supplies a unified description of 
direct and resonance reactions. 

Reaction channels that are energetically allowed are called open channels, and all 
others are closed channels. In the open channel /3m, the wave function at large sep- 
arations of the /3 fragments assumes the form 

< 10 ) 

A projection P which selects the open channels is a projection operator for which P\l> 
has the same asymptotic form in any of the open channels as ^ itself. Thus, at large 
particle separations in the open channel /3m, P\I> must yield 

P ^ f /3m^/3m ( 20 ) 

The open-channel subspace into which P projects will be called P. Since only the as- 
ymptotic behavior of P# has been uniquely prescribed, the subspace & is not unique 
and, hence, neither is the projector P. The largest subspace SI orthogonal to & can 
be referred to as the closed- channel subspace, and its projector is Q = 1 - P. Since the 
incident channel and all other open channels are included in P, when Q operates on a 
state S&, it yields a function which in any channel has the asymptotic behavior 

Q^/hn 0 for any i3m (21) 

In configuration space, Q'J' is limited essentially to the regions in which fragments of an 
open channel are interacting. 

The pair of coupled equations obtained by operating on equation (5) in turn with P 
and Q is 
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(e - jr pp )p* = ^pqQ^ 


(22a) 


(E - ^qq)Q^ = -*q p P* (22b) 

Li view of the boundary condition (eq. (21)) that Q 1 !' satisfies, the formal solution of 
equation (22b) can be written as 


Qtf = 


E - #>, 


QQ 


•*Qp P * 


(23) 


When equation (23) is substituted into equation (22a), a Schrbdinger equation for P^ is 
obtained 


(E - H)P# = 0 (24) 

with the effective Hamiltonian H given by 


H - ^p p + **pQ — -*Qp (25) 

~ ^QQ 

The transition amplitudes obtained from P$r are exactly those given by the full state 
vector i I r ; thus, resonance effects occurring in the reactions are accompanied by rapid 
energy variations in the effective Hamiltonian H that governs the behavior of P#. It is, 
in fact, the second term on the right side of equation (25) from which the rapidly varying 
part of H arises. The compound states of the total system are described approximately 
by the eigenstates (p v of ^qq belonging to the discrete spectrum 

( g v - = 0 (26) 


These eigenstates are true bound states and, hence, have an infinite lifetime. In addition, 
there is the continuous spectrum of eigenvalues that begins at a minimum value equal to 
the smallest threshold energy for opening any of the channels which are closed at the en- 
ergy E and which are, thus, not in the set of open channels selected by P. Solutions of 
the homogeneous form of equation (22b) are 6(E - but any admixture of these 

functions to the solution given by equation (23) is ruled out by the physical requirement 
that 
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(<P„|<») = (E - V 1 <V„|j«’qp|p*) 


be well behaved as E — S v . 

The exact compound state is described by OF and decays into the open channels by 
virtue of the coupling provided by the term on the right of equation (22b). But, for values 
of E in the neighborhood of an isolated resonance, the compound state has a long life- 
time and the rapidly varying part of Q’k is proportional to the state <p p whose eigen- 
value £ is close to E. At an energy where resonances overlap, a linear combination 
of the states cp v , whose eigenvalues lie within the region of overlap, approximates the 
compound state. These properties are verified by expanding the operator (E - -^qq) ^ 
in the total set of eigenstates of including those in the continuous spectrum, 

and by then using the expansion in equation (23). When the same form of the operator is 
employed in equation (25), H can be separated at an energy value E, where resonance oc- 
curs into two parts, one of which, H', is slowly varying, and the other, Hj^, is rapidly 
varying with energy: 


H = H’ + H r (27) 

The resonant part Hj^ consists of a single term from the expansion when E is near an 
isolated resonance. When resonances overlap at E, consists of a finite sum of 
terms corresponding to the contributions of the states cp v whose eigenvalues & v lie in 
the region of overlap. The complete derivation of equation (27) is given by Feshbach 
(ref. 2). The main elements of the development are sketched in appendix B. 


Appropriate Forms for Projector That Selects All Open Channels 


Consider now an open channel am. The projection of 'k into the subspace 

However, the full 


am 

T'" certainl y approaches \l> itself asymptotically in channel am. 


projection into is not required to obtain the same asymptotic behavior; that is, 

there exists, in general, a subspace 'f'iZl of such that the projection of ^ into 

this subspace still yields the correct asymptotic result in channel am. Let ' 

designate the largest subspace of . which is orthogonal to Then the projec- 

f /Q wIH ^ HI 

tion of f in f 


(/3,C) 
am 


vanishes asymptotically in channel am: 


C) 


am am am 


(28) 
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The corresponding projection operators are 


n = 

am am 


+ n 


(/3,C) 

am 


where 


n ( a ) JjO 3 ? C) Q 


am am 


Thus, a decomposition of the subspace '^ am is introduced which is not necessarily 
unique but which has the following properties: 


am am 


(29a) 


n ^’ c ^— >0 

am am 


(29b) 


and 


n am«« = 0 


(29c) 


In order to see that ^ is generally a somewhat larger projective space than 
necessary in which to retain the desired asymptotic information about in channel am, 
other possible requirements can be cited which the subspace ^ may have to fulfill. 
First, consider an open rearrangement channel j3n in which some fragments of the set 
a are bound together and the others are the same, as for example in neutron capture. 
When the internal state cp^ n is expanded in the set <P am , clearly, large components for 
open channel values of am can occur. But the contribution of '^ am that enters into 
the description of channel /3n corresponds to functions that vanish asymptotically in 
channel am. The part of ‘f ' ^ that yields only functions of this kind cannot enter into 

the asymptotic behavior of 'k in channel am. However, it contributes to that part of & 
which plays a role only in the description of rearrangement channels. 

A second function that might conceivably be required of involves its projec- 

tion in 2. A restriction that is always implied on the extent of S? is that the closed 
channel subspace Q must have enough in it to provide a description of compound states. 
The compound states in nuclear reactions, however, are expected to have a very small 
overlap with the relatively low lying states <p available to inelastic excitation of the 
target system. In the case of an incident nucleon, the compound nucleus has an excitation 
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above its ground state of the order of 8 MeV plus the kinetic energy of relative motion. 
The actual complex motion of the compound system can be pictured as reached through a 
succession of two-body interactions. It has been estimated (ref. 14) that, even at low in- 
cident energies, this multiple excitation produces about seven quasi- particles (i. e. , eight 
excited particles and seven holes compared with the initial state of the target). Thus, the 
portion of Y am required for this second purpose is very small but not identically a null 

space. In both the examples presented, the parts of Y„ m that were utilized meet the 

(B Cl “ ni 

criterion established for . 

To illustrate how the separation of Y am might be effected in terms of a set of basis 
vectors in the subspace, consider a potential function in the mass center coordinates of 
the a fragments with the same symmetries in these coordinates as Jf itself. For nu- 
cleon scattering, a shell model potential for the target nucleus would be appropriate. In 
equation (8), the coefficient f can be expanded in the complete set of eigenfunctions 
of the potential with suitable outgoing wave boundary condition. For large separations of 
the mass centers, the discrete bound eigenfunctions do not contribute to the expan- 
sion. Thus, can be taken to be a space spanned by all or some subset of the 


basis vectors 

The subspaces Y am , and '^'am > are ort hogonal to all other Y an , where 


,(a) 

am’ 


4P, c) 

am 


n * m. Unlike Y am , however, the orthogonal subspaces Y^^ and 


are not in- 


variant under the operation of K^, which does not commute with either II or 
Let the set of all open a channels be designated as A. These open channels are enu- 
merated by the indices am, in order of increasing values of the internal energy jt a , 
when m runs from 0 to N(a). Also, let the vector sum of the subspaces Y am over 
the set A be called Y^, the subspace orthogonal to it be Y^, and the projection oper- 
ator into Y^ be n^. Then these quantities can be written in the form 


*A = ^al® ^2® * 


® y 


aN 


•5^^ _ /T\ ^ £T\ 

' A ~ a, N+l® a, N+2® * 


(30) 

(31) 


and 


A 

^A = ^ v "am 
m=0 


(32) 


and similar expressions for the vector sum Y^ of the subspaces Y^^ over the set 
A, and its projector etc. can be written. 
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Consider the intersection of 9 and Y^, which is a linear subspace of both 9 and 
For an arbitrary projection P, this plane of intersection will be skew to the ortho- 
normal basis for Y^ constructed from the eigenstates of Jt a (each vector of the basis 
set has the form k<P am where m c A). However, there is sufficient flexibility in the se- 
lection of a projection into the open channels that in any problem P can always be chosen 
so that 


9 n y a = Y^ 


This follows because the part of 9 required for the description of the asymptotic be- 
havior in the open a channels must also lie in Yp^ and, by definition, the subspace 
is adequate for this purpose. The subspace in which is orthogonal to 


Y^\ is the part of 9 which is required only for the description of open rearrangement 
channels. (The space Y^ ' may also contribute to rearrangement channels involving 
bound a fragments.) These subspaces are visualized schematically in figure 1. Thus, 
P can always be chosen to have the form 




p = n 


n( a )n 
ll A u o 


03 ) _ 


= 0 


(33a) 


with the projective space 


9= Y^@Y^ (33b) 

Such a choice has the advantage that the projective space Y ^ of 11^, which yields 
the projection into the open a channels, is spanned by an orthonormal basis whose vec- 
tors are all eigenstates of jt a . This choice is convenient for separating the elastic and 
inelastic processes from the remaining reactions. However, none of the spaces in equa- 
tion (33b) is invariant under K a and none of the projectors in equation (33a) commutes 
with K a . 

From the discussion given previously, the assumption that the compound states have 
no components in Yp^ does not appear to be overly restrictive in nuclear reactions. It 
follows from this assumption that none of the subspace Yp^ is needed in 21 and, hence, 
that 9 can be enlarged to include all of Yp^. Then the form 

p = n A + n Q , n A n 0 = o (34a) 

can be chosen with the corresponding relation 
14 



p = 


Yp@Y 


O 


where Yq is the subspace in & which is orthogonal to Yp^. Although IIq in equa- 
tion (34a) does not commute with K^, 11^ (eq. (7d)) does. 


Two Kinds of Open Channels 


In order to make more explicit the general arguments given so far, a somewhat sim- 
plified reaction where only two kinds of channels are open, the a and the /3 channels, 
is considered. Let A and B designate, respectively, the complete set of open a and 
open /3 channels, and let and Ilg be the projectors of the sort defined by equa- 
tion (32) into the respective total subspaces Yp^ and Y^, like those shown in equa- 
tion (30). Then, with the assumption that neither Yp^ nor Y-q contribute to Q, under 
which projectors P of the type given by equation (34a) are valid, P can be chosen to be 
the projection into the space *A+B spanned by all the vectors of Yp^ and Y^ com- 
bined; that is, by defining the subspace 


*A + B = 



(35a) 


the projector P is chosen to be: 


P = n 


A+B 


(35b) 


An alternative prescription to that of equation (35a) for the subspace ^ + b can be 
obtained from the following constructive argument. Resolve all vectors of Y-q into two 
components, one contained in Y^ and the other orthogonal to Y^. The set of compo- 
nents in Yp^ span a subspace (lying entirely in Yp), which is the projection of Y^ into 
Yp^. The set of orthogonal components spans a subspace, called Y-q^, that is orthog- 
onal to Yp^. Then is the vector sum of Yp^ and Y-q because any linear 

combination of a vector in Yp^ and a vector in Yq can be expressed as the sum of two 
orthogonal vectors, one in Yp^ and the other in Yq p^: 


^A+B = *A® *BpA 

= ^B® %B 

(36) 

n A+B = n A + n BpA’ 

n A n BpA = 0 

(37) 

= n B + n ApB’ 

n B n ApB = 0 

(38) 
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If the incident channel is of the a type, the form given by equation (37) corresponds to 
equation (34a), but, if it is a /3 type, equation (38) corresponds to equation (34a). It de- 
pends on the details of a given reaction whether or not Il^'t arising from the use of 
equation (37) contributes to the asymptotic behavior in the fi channels, and likewise for 
Ilg’k in the a channels when the form of equation (38) is used. 


Actually to construct projection operators like n 


BpA 


in a given problem may not be 


easy, as can be appreciated from the example of stripping or pickup reactions discussed 
in appendix C. A possible simplification in its use is seen upon writing 


n BpA^ " II BpA II B' ir + n BpA^ " n B^ 


and noting that only the first term on the right enters into the asymptotic contribution of 
n BpA * in the open /3 channels. The projected vector (1 - IIg)'l> contains no such as- 
ymptotic information whatever, nor does any further projection of this vector. Equa- 
tion (37) shows that 


n A+B n B " n B n A+B " n B 
and, after inserting equation (37) for n^ +B , 

n BpA n B * » - n A> n B < 39 > 

is obtained. This operator is not itself a projection operator, nor does it yield one when 
added to 11^. (The only exception to this statement occurs when and Ilg com- 
mute, which corresponds to the special case that is orthogonal to and, hence, 

is the same as ) The operator (eq. (39)) is easily constructed, however, and 

when solving equation (24) only that part of llgp^’k given by it needs to be obtained. 

When more than two kinds of channels are open, the form of P given by equa- 
tions (37) and (38) can be generalized in terms of projections into additional orthogonal 
subspaces. The alternative expressions for P reveal that a projector of the form of 
equation (34a) can also be written as 

p = n B + n Q „ n B n G , = o 
= n c + iIq,,, n c n 0 „ = o, etc. 

where n B is the projection into the open j3 channels, II ^ is the projection into the 


I 
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open y channels, etc. These alternative forms of P are not required here. They cor- 
respond to a multicomponent form of the projector given recently by Hahn (ref. 15) for 
application to multichannel processes. 


Generalized Optical Potentials 

The generalized potentials that describe the scattering into a selected set I of a 
channels, which includes the incident channel ao as well as an arbitrary group of in- 
elastic channels, may now be derived. For the projection operator P into all the open 
channels the one given in equation (34a) is used. The projector n A defined by equa- 
tion (32) can be resolved as follows: 


n A - n i + n A I 

(40a) 

n i n A-I = n i n O = n A-i n O = 0 

(40b) 


Both IIj and n A _j commute with K fl! but as already indicated the projector IIq does 
not. Combining equations (34a) and (40a) yields 

p = n j + iijr, rijH^ = o (4ia) 


where 


n K * n A-I + n O < 41b > 

is the projector which selects all the open channels except those in the set I. The K 
channels comprise the excluded open inelastic channels as well as the open rearrange- 
ment channels. 

In Feshbach's theory, equation (27) expresses the effective Hamiltonian H for P'l' 
as a part H’, which is almost energy independent, plus a resonance term H^. Further 
details are reviewed in appendix B, where H’ is given by equation (B2) and a suitable 
expression for Hj^ is developed. In terms of the incident channel resolution of # 
given by equation (1), H’ has the form 


H* = PK P 

a 


+d p 


where 
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Up = PV a P 


+ U T 


and U D represents the effect of all distant resonances. 

The generalized potential or effective Hamiltonian for in a form which separates 
the slow from the rapid energy dependence may now be obtained. For this purpose, 
equation (24) is written as 


(E - H')P^ = H R P'$ r 


(42) 


Operating on equation (42) with II K and solving for 4> K s yields 


*K = 


E " H KK + 


< H fcl*I * ^P*) 


(43) 


where 


H ki= n K H,II p etc - 


The first term on the right arises from the coupling between the I and K channels, and 
the second represents the direct coupling from the resonance states into the K channels. 
When equation (43) is used with the projection of equation (42) by IIj, the result can be 
written as 

(e - h^ (44) 


where 


M) = 


= H n 


H k 


E - 


h kk± 


H kl 


(45) 


and 




(±)_ 


= n x + n 


K 


E - 


H ki 


(46) 


so that 


h^ = 
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muni hub i« 


In order to interpret equation (44), the analogous operations are performed on the equa- 
tion 


(E - H’)*^ = 0 

where ^rn) is t ^ ie ou tg°i n g wave solution with incident wave in channel a Q (am), 

with the result that 


^oK s n ^ (±) = 


JTo 


E-H^ + fc 


H* i i/W 
h KP'oI 


(47) 


and 



Thus, 


• n^Sf’ 

is the eigenstate of h^ with the same asymptotic behavior in the I channels as 
itself. Also, equation (47) shows that 

^ +^ ( 0 ^ = n (± V$ (49) 

The states i (^ia) are or ^honormal, but their projections ('/'ml) are no ^' 

Equation (44) for is in the desired form. The resonance contributions to 
are isolated on the right side. The first term on the right gives the direct contribution 
into the I channels from the resonance states, and the second term corresponds to the 
indirect contribution that arises as a result of the resonant behavior in the K channels 
and their coupling to the I channels. Equations (18) and (44) are completely equivalent 
to one another, but equation (44) has the advantage of exhibiting the resonance contribu- 
tions explicitly. The expression for obtained by solving equation (44) is in the form 
of a slowly varying part plus a rapidly varying part: 

*1 = ^ + - (50) 

E - h( + ) + ie 
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The transition matrix, which follows from equation (50), is 


T(I) = T p (I) 




+ ( 4 ' 


,(-) 

ml 


O^'^Hb I P* 


R 


*> 


( 51 ) 


The transition matrix element T p (I) given by xp^ is clearly equal to the matrix ele- 
ment Tp given by xp^ for any of the I channels, because xp^ has the same asymp- 
totic behavior as xp^ in these channels. The total transition matrix element (eq. (50)) 
is the same as that which follows from P 1 ^ (eq. (B4)). This correspondence can be veri- 
fied immediately with the help of equation (49), which gives 



= n 



(-) 

ml 


In order to obtain the generalized optical potential for inelastic scattering, the tran- 
sition amplitude (eq. (51)) is averaged over a suitable energy interval AE, as described 
in detail in appendix B. The steps involved are exactly analogous to those which lead 
from equation (B4) to equation (BIO), and the result is 


T(I) = T p (I) 


iir(xp 


i77 ^ml 


•^pqA 


(+) 




(+)' 

ol 


(52) 


where 


A 


A 

AE 


and 


A = 



and where the sum runs over all discrete eigenstates <p ^ whose eigenvalues 
within the interval AE. 

The optical model Hamiltonian 



lie 


= h(+) 


h M= h 


+ U, 


CN 


(53) 
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can be characterized as follows: At energy value E, the eigenstate 0 ^ + \ which satisfies 
the same boundary conditions as 'i'p gives the transition matrix (eq. (52)). An integral 
equation satisifed by eigenstates of h M is 


0 (+) = *<+> 


- h^ + ) 


U CN 0 


(+) 


E - h 


+ ie 


while a state with the desired transition matrix is 

« <+) - - *” 7 t — (a ( -> VpqAjfQpQ*^)^^ (54) 

E - h' + ' + fe 


Comparing these last two equations shows that 


U CN ~ _ur ^ 

1 


hr 


1 

E - h^ + ie 


1 


(55) 


An alternative approach to the desired optical potential is to obtain it directly from 
the optical potential H M whose eigenstate yields the energy- averaged transition 
amplitudes for all open channels: 


(E - H m )x (+) = 0 (Bll) 

The potential has been derived by Feshbach (ref. 2) and is shown in equation (Bll). 

From equations (B12) and (B15), the generalized optical model Hamiltonian for all chan- 
nels can be written as 


H M = PK a P 


+ U 


where the generlized optical potential U is 

U = p v a p + U D + U CN 


and 
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1 


u, 


CN 


17T 

AE 


J^PqA 


A*. 


1 - — AWA 
AE 


QP 


The last equation is a more symmetrical version of equation (B15), and W is defined by 
equation (B7b). 

It is easy to project out of x^ + the part that describes the scattering in the selected 
a channels and to determine the optical model potential which governs this part. It fol- 
lows as a consequence of there being no closed channels in the optical model description 
that 


PX (+) = X (+) 


and that 


ph m - h m p - h m 

as is easily verified from the explicit expressions for x^ and H M . The form of P 
given in equation (41a) leads to the resolution 

x< + > = 4 + > + xk +) 


Operating on equation (Bll) with II K and solving for Xj^ give 

1 


v(+)„ 

X K 


E-H^ + ie 


H M x (+) 

n KI X I 


(56) 


where = n„H, T II T , etc. Using equation (56) in the equation obtained by operating 

Jtv JVL I 

with IIj on equation (Bll) then yields 


|E “ Hji " < 


1 


E - + * 


Elx< +) = 0 


H ia) x i 


(57) 


which is the desired expression. 

Two expressions (given by eqs. (53) and (57)) for a generalized optical potential des- 
cribing the scattering in a selected set of inelastic channels have been derived. The two 
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different forms were obtained by reversing the order of the operations of energy averag- 
ing the transition amplitude and projecting into the selected channels. There is no 
a priori reason for these to be identical potentials, since the condition imposed upon 
them is only that their corresponding scattering states have the same asymptotic behav- 
ior. However, they do, in fact, give the same potential. Their equivalence is most 
easily established by showing that their scattering states are the same everywhere. The 
proof of this follows directly by operating with the projector nj on equation (B13b) and 
making use of the identity 


n. I 1 o<-)t 

+ k E-h< + > + t£ 

Comparison with equation (54) then shows that xj + ^ = 0^. 

In the development of this section, the commutativity of IIj and has not been 
used as yet. Hence, the equations have the same form when equation (33a) is used for 
the projector P. The use of equation (41a), however, does allow for additional simplifi- 
cations. For example, the terms in equation (57) can be written as 


H n - K « n i + u n 


TjM _ T T Tiki _ II 

H nc ~ U IK’ H KI “ U KI 


H KK “ n K K o; n K + U KK 


with Ujj^ = njUn K , etc. Equation (57) then becomes 

1 


E • K a - D n ' D nc 


E - H KK + i£ 


u kiM +) = o 


PROPERTIES OF GENERALIZED POTENTIALS 

Because PSk is the projection into all open channels, its effective Hamiltonian H 
is Hermitian as are the separate terms H’ and H^. Inspection of equations (25) and 
(Bla) immediately verifies that these operators are indeed Hermitian. The slowly vary- 
ing part H’ contains in the terms PV^P and Up the interaction mechanism that 
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causes direct transitions to the open channels. The term Up represents the net effect 
of all virtual transitions from the open channels to the closed and back. The generalized 
potential in equation (44), which governs v&j, is not Hermitian, however, as follows from 
the imaginary term ie in the denominators of equations (45) and (46). The imaginary 
terms represent an absorption of the incident particle arising from processes in which 
the exit channel belongs to the set of excluded open channels K. 

The physical significance of the individual terms appearing in equation (44) may now 
be interpreted. The part h^, which is slowly varying with energy, is the contribution 
to the potential from the direct processes. Equation (45) can be written in the form 

h< +) = HJj ♦ Pr -X_ Hja - - H^HJg 

" h kk 

where Pr stands for the principal value. Because HJj and are Hermitian, the 

first two terms are Hermitian and the third is imaginary. The first term contains 
njV^nj and, in addition, ITjUpIIp which is the result of virtual transitions from the 
open I channels to the closed channels and back to the I channels. The second term is 
the off- energy- shell contribution of the propagator and corresponds to virtual transitions 
through the closed channels, which are connected to the I channels by way of intermedi- 
ate states belonging to the open channels K. The third term comes from the energy- 
conserving part of the propagator and represents the effect of the direct transitions to 
the set K of excluded open channels. 

The resonance terms containing H^ arise from the reaction mechanism which pro- 
ceeds via real transitions to the discrete compound states of the closed channel subspace. 
The term IIjHjjP^ gives the decay from the compound states directly into the I chan- 
nels. The imaginary term from the on- energy- shell part of the propagator in equa- 
tion (46) involves the decays into the K channels. In both the direct and the resonance 
terms, the imaginary part will tend to increase when the number of inelastic channels 
assigned to the set K is increased, because the transitions associated with the propa- 
gator term will become more numerous thereby. 

The generalized potential v< + > has many of the properties of the single- channel po- 
tential derived by Feshbach (ref. 1). Only slight modifications of his discussion are re- 
quired here; therefore the discussion is limited mainly to stating the corresponding prop- 
erties. The smallest threshold energy for opening any of the channels not belonging to 
the set I is e^-. The abbreviations -^qq - OjJfOj, = IIj-V^Oj, etc. are also 
adopted. The eigenstates of -^qq are displayed as follows 

» 

jf ©©V = e m T ?m 


24 



for the discrete bound eigenstates and 


•*e 0 n(e y) = e r) 


for the continuum eigenstates. The continuous spectrum begins at the minimum value 
e K . The additional labels y are those used to designate different states with the same 
eigenvalue e'. 

The generalized potentials H and are nonlocal. This property can be exhib- 

ited explicitly by expanding the inverse operator of equation (25) in the complete set of 
eigenstates of -*qq and the inverse operator of equation (17a) in the eigenstates of 
In the case of Vj + ^ the result is 


Vj +) (E) = Vg + 


I 


T 01 

10 


E - e„ 


_ r Q ! 

v 0i + 


m 


f*f 

J Je K 


de’ V- 


a 
1 9 


E - e ’ + fe 


t ol 

7 ei 


(58) 


In treating the energy dependence of Vj , the set of channels I is kept fixed, even at 
small values of E where some of the channels in I are closed. The derivation and 
formal expression for y( + ) 

are the same at any value of E. However, it should be 
noted that the interpretation of equation (44), as the separation of Vj + ^ into a slowly 
varying term and a resonant term, becomes invalid when IIj includes closed inelastic 
channels whose projection subspace enters significantly into the description of the com- 
pound states with energies close to the value E. 

From equation (58), 


ReVj + ^(E) = 


tOl hir^^m 


10 


E - € 


^I + 


m 




+ Pr 


f f ^ v i© *» y))(^ ’> y) | v oi ( 59a ) 

Je K E - e 'J 


and 
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ImV< + >(E) = -7T f dyV^ |t 7(E, y)) ( V (E, y)| vgj if E 


> € 


K 


(59b) 


= 0 if E <e 


K 


From these expressions, it is clear that the imaginary part of Vj + ^ is nonpositive, as it 
must be if it is to represent an absorption, and that the energy derivative of the real part 
is nonpositive. A dispersion relation between the real and imaginary parts follows di- 
rectly from equations (59a) and (59b): 



The optical model Hamiltonians given by equations (53) and (B12) are, of course, 
not Hermitian. The spread in energy introduced by the energy averaging allows for an 
interpretation in terms of a time-dependent wave packet picture. The optical potential 
describes only the direct or prompt processes accompanying the scattering of the wave 
packet. The time-delayed part of the wave packet resulting from the formation and sub- 
sequent decay of the compound system is not contained in the solution of equation (Bll). 
The averaged effect of all the resonances on the direct reactions is given by Uq N , the 
imaginary part of which accounts for the attenuation of the packet due to the effective 
absorption of direct particle flux into the compound- state reaction mechanism. In the 
limit that the matrix elements of W are small, is pure imaginary. 

The term u CN given by equation (55) similarly serves to give the effect of the res- 
onances in the generalized optical potential (eq. (53)). Its form differs from that of 
because it operates on the projected wave packet, which describes the scattering in the 
I channels only. However, the compound states are reached by transitions emanating 
from both the projected part and the projected- out part of the wave packet, and, likewise, 
the decay of the compound states feeds into both sets of channels. The form of u^ N 
takes account of these effects. It is also clear that H M and h M are nonlocal potentials. 
In the set of coupled equations for the relative motion in the inelastic channels, the ma- 
trix elements 
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enter as coefficients. These are complex functions of the relative coordinates. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, February 14, 1967, 
129-02-07-07-22. 



APPENDIX A 


A 

d 



c n< r o> 


D 

E 

AE 

E i 

f i< F d> 

f /3m 

G(E) 

g(r d , V) 

g(r 0 ,fi) 


H 


H’ 

M 

U JK 

H k 


H 

H 


M 


R 


SYMBOLS 


set of all open a channels 

target nucleus consisting of nucleons 1, 2, . . . , A 
set of all open /3 channels 



composite particle consisting of nucleon 1 bound to nucleon n 

total energy 

energy interval 

eigenvalue of m (*^qq + W)M 

discrete eigenvalue of ^qq 

V 

coefficient of <P am in expansion of ^ in eigenstates of 
coefficient of <Pp m in expansion ^ in eigenstates of 
Green's operator 


w-oKy r 

effective Hamiltonian for 

part of H which varies slowly with E 


n i H ' n K 

optical model Hamiltonian for all open channels 
part of H which varies rapidly with E 
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total Hamiltonian 

•*PQ 

P.JTQ 

*^©0 

epjfBj 

h (±) 

O^HTIt 


hj(? 0 ) (b|Yj* d ) 

optical model Hamiltonian for selected inelastic channels 
I selected set of open a channels including incident channel ao 

K set of all open channels except those included in set I 

K a noninteracting Hamiltonian for channel a 

X a sum of internal Hamiltonians of particles in a channels 

k specifies wave vectors of noninteracting channel fragments 

k Q specifies wave vectors of incident particles in incident channel ao 

k m specifies wave vectors of scattered particles in exit channel am 

M(E) projector into subspace i2^(E) 

Mj mass of particle i 

M normalization constant for vector z" _ 

mn mn 

N number of resonances in energy interval AE 

N normalization constant for vector 4>’ 0 

N(E) dimension of subspace ^Jp(E) 

N(a) number of open a channels 
n designates internal states of channel fragments 

n nucleon 

P projector which selects all open channels 

P F n FB + n ApFB 

& projection subspace of P 

Q 1 - P, projector into closed- channel subspace 

Q, projection subspace of Q 
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R/ 


D 


r i 


(E) 


T(I) 

T(amk m |aolT o ) 


a 


U 

U 


CN 


U, 


u. 

U- 


D 

(+) 

I 

IK 


U 


u 


CN 


u 


Q!k 


vJ +) 


projection subspace of M(E); space spanned by all 4> v whose levels 
£ overlap at energy E 

an effective radius for nucleus A 

position vector of particle D relative to center- of-mass of total system 

position vector of particle i relative to center- of-mass of total system 

transition matrix element for scattering into set of channels I 

transition matrix element connecting incident channel ao to exit chan- 
nel am 

sum of center- of- mass kinetic energies of particles in a channels 
transition matrix element which follows from H' 
generalized optical potential 

contribution to from averaging over resonances in energy inter- 

nal AE 

contribution of distant resonances to H’ 
operator satisfying V a ^f = 

n^K 
pv^p + u D 

contribution to h ffi from averaging over resonances in energy inter- 
nal AE 

eigenstate of T a which is product of plane waves with wave vectors 
specified by k 

generalized potential for selected inelastic channels 


V K> 

n i V a 0 I 



V a 

interaction between particles in a channels 




1® 2® * ’ * vec tor sum of subspaces 

r over set A 


*a,N+l®*ot,N+2® ■ • • sum of subspaces 

in r A 

r 

am 

not contained in 

*1°° 

r al sr a2 B ' ■ ■ 



*A + B 

subspace spanned by all vectors of and 


combined 
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y'O 3 , c) 
am 

W 

X t (r D ,rf) 

X i( ? o) 

x i 

x ! 

x r 

Y j< F D> 

y j 

y! 

y l 

Z mn( ? D^> 

z mn 

Z M 

mn 

a,/3 

am 

/3m 


largest subspace of *B which is orthogonal to Y ^ 

space spanned by all vectors of the form Z(rj^,^)i^ o 

direct sum of Y m inset I 
am’ 

subspace in ^ orthogonal to Y^ 

subspace in & orthogonal to 

subspace spanned by all vectors of form X<P am 

subspace of Y am such that yields asymptotic behavior of ^ in 

channel am 

largest subspace of orthogonal to 

■X' Q p(E-H' + i£)-V pQ 

X i< ? 0> 

complete orthonormal set of functions of r o 

X i<*>o 

(l - n B ) Xi 

N(1 ^FB^ x i 

complete orthonormal set of functions of r^ 

Y j( r D> 

Wo 

(1 - n A > yj 

complete orthonormal set of functions of r^ and rf 

Z mn<?D’V) 

Z mn^ o 

M mn^ " n A^ z mn 

indices specifying types of channels 

indices specifying the various a channels 

indices specifying an exit channel of the /3 type 

- ImE. 

2 1 
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€ 

e K 

V 

v(e , y) 

V 

* 1^0 


0 

e 

A 


I 

(+) 


A 

V«i> 


positive infinitesimal 

smallest threshold energy for opening up any of the channels not in the set I 
position vector of particle i relative to center of mass of nucleus si - 1 

continuum eigenstate of -#q q 
discrete bound eigenstate of 

position vectors of particle 1 and incident nucleon respectively from center of 
mass of nucleus D 

1 - n T 


continuum eigenstate of h. 


m 


E 


AE 
(AE) _1 A 

complete orthonormal set in variable 

position vector of particle i relative to center of mass of nucleus si 


n A 

projection operator into subspace 

*A 

n (°0 

n A 

projection operator into subspace 

r A 

n A-I 

n A-"l 


n A+B 

projection operator into subspace 

*A + B 

n BpA 

projection operator into subspace 

^BpA 

n FB 

projection operator into subspace 

^FB 

n l 

y v n «m 



me I 


n K 

n A-I + II o 


n o 

projection operator into subspace 

*o 

n w 

projection operator into subspace 

r (P) 

G 

n am 

projection operator into subspace 

am 

II^) 

am 

projection operator into subspace 

iA a ) 

am 
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n cfm^ projection operator into subspace 


■WC) 


n 

eigenstate of jf 



eigenstate of 


*o 

ground state of nucleus 

J 

*o 

N &o " g( ? l)^o] 



discrete eigenstate of 

W 

X 

any state of motion for mass 

X (+) 

continuum eigenstate of 

h m 

*(+) 

X I 

n r x (+) 


X W 

X K 

n y( + ) 
11 K X 


* 

Q!OK 



©jtf 


*1 

iij* 



*<+V 

Q!Ok 

fa) 

*i 


(±) 

m 

o 

w 






CO 


(±) 


continuum eigenstate of 
ground state of nucleus D 

discrete eigenfunctions of potential function of mass- of- center coordinates of 
a particles 

complete set of states of nucleus D 

continuum eigenstate of H* 
ground state of nucleus sd - 1 

continuum eigenstate of H’ 

V? 

n i+ n K „ J H ki 


E " h kk * ie 


M(E - *QQ>’ ljr < 


Qp 




f ) 
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APPENDIX B 


REVIEW OF PERTINENT CONCEPTS FROM FESHBACH'S 
UNIFIED THEORY OF NUCLEAR REACTIONS 

When the inverse operator in equation (25) is expanded in the complete set of eigen- 
states of -*qq, the energy dependence of H in the neighborhood of E can be expressed 
as the sum of a slowly varying part and a rapidly varying part: 



M 


In the first term on the right, which is called Hj^, the sum extends over the discrete 
levels <? , which are sufficiently close to the value E that their contributions vary rap- 
idly. The term Uj^ represents the remaining terms from the sum over the complete 
set. It is the net contribution of all virtual transitions from the open channels to the dis- 
tant resonances and back and varies slowly with E. The strongly energy- dependent part 
of H can be written as 

Hr = ^pqM(E) _ 1 — *- Qp (Bla) 

* ' "^QQ 

where 

M < E >= Z | *„)<*„ | < Blb > 

The set of p values in the summation, called N(E) in number, depends on the value 
of E. In the region of isolated resonances, N(E) is at most equal to 1. Note that M(E) 
is the projector into the N(E) dimensional subspace ^j^(E) of SI spanned by the corre- 
sponding set of orthonormal vectors <p The remaining part of H is the slowly vary- 
ing part H’, so that 

H' = JTpp + U D (B2) 
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The integral equation for P# in terms of these two parts of H is 


FF = ,//(+) + 1 H n P$ (B3) 

° E - H» + ie K 

where is the eigenstate of H' with the same boundary conditions as 1 3>. The am- 

plitude for the transition from the incident channel ao to the open channel am is 


T(amk m | aok 0 ) = T p + (^ } H R | P*) 


(B4) 


Here, T is the transition matrix element 


where k Q and k m designate, respectively, the incident and final wave vectors of the 
a particles and is the eigenstate of H’ with the energy value E, and the boundary 

conditions of incoming spherical waves in the open channels and plane waves in channel 
arm with wave vectors specified by k m , 
which follows from i ]/^ . 

There is a rapid energy dependence from the second term on the right of equa- 
tion (B4), and it originates from both H R and P'k. (It is assumed that the energy is not 
in the neighborhood of a threshold for opening a new channel, so the projection operator 
P does not change with energy. ) In order to express the resonance variation explicitly, 
the rapidly varying part is written as 


H r P* = 


■*PQ 


u> 


and a more suitable form of the vector 

co = M(E)(E - jTqq) “ 1 ^Q p | P*) (B5) 

is obtained. From equation (23), co is seen to be the projection of Q’F into ^ R (E). 
When equation (B3) is used in the defining equation (B5) and the commutativity of -^qq 
and M(E) is used, an equation for o> results that can be expressed as 

(E - ^q Q - MWM)w = M -*Q P ^ +) (B6a) 


where 


W = JT, 


QP 




E - H» + ie 


PQ 


(B6b) 
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Equation (B6a) is a linear equation in the subspace ^j^(E) and is conveniently solved in 
terms of its Green's operator: 




E ,*E, 


(B7a) 


Eh)K|- EkX'v 

i IX 


(B7b) 


Then the desired expression for is 

co = G(E)^ Qp ^ +) 

for, when this is inserted into equation (B4), the transition amplitude becomes 



i 


In this form, the energy variation of T in the neighborhood of E arises almost 
entirely from the demoninators in the second term. Let r. be one-half the imaginary 
part of Ej. It is assumed that I\«Re(E.). It is also assumed that E. and co. have 
only a slight energy dependence in an interval of the order of Tj about the energy value 
Re(Ej). Then it is meaningful to speak of an individual resonance E., even in the over- 
lapping resonances region, and the energy width about the value Re(E-), in which it 
contributes a rapid variation to the transition amplitude. 

The transition amplitude is now averaged over an energy interval AE which is big 
enough to encompass a large number N of resonances and to be much greater than r., 
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but yet is sufficiently small that T p (or equivalently H’) is essentially constant over the 
interval. (The term U D in H' varies slowly with energy over an energy interval con- 
siderably larger than the average value , because only the "nearest" distant reso- 
nances change with E and these constitute a small potion of all the distant resonances 
contributing at any given energy. ) It follows from the assumptions made that to a good 
approximation 



(B9) 


The sum over resonances now includes all lying in AE. As E varies over the 
range of integration, a cumulative vector sum of the projective spaces ^? R (E) gives a 
total projective space spanned by all the eigenstates cp^ with values of laying in 
AE. The projector for this subspace associated with AE is 

A - EkX’vl 

AE 

Under the assumption that the uk are energy independent, both sums in equation (B8b) 
may be extended to run over AE which permits equation (B9) to be replaced by 

t ■= T P - («&' | w*qp K +) > <bio > 

It is clear that the form of expressions like (BIO) is the same for either the case of iso- 
lated or of overlapping resonances. 

The generalized optical model Hamiltonian H M for all the open channels is defined 
by the condition that at the energy E equal to the average in AE, its scattering state 
X (+) , which satisfies 


(E - H m )x (+) = 0 (Bll) 

with the same boundary conditions as P'k, yields the transition matrix T. The part 
Tp of T corresponds to the scattering given by a state of H', so that an additional part 
Ucn i s appended to H', in order to produce the contribution arising from the averaged 
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effect of the resonances, that is, the second term on the right of equation (BIO) 


H M = H ' 


+ u, 


CN 


(B12) 


An integral equation satisfied by is obtained from the two-potential form (eq. (B12)) 

of H m , namely, 


X (+) = 


E - H’ + ie 

,(+> 


U CN X 


(+) 


(B13a) 


On the other hand, an explicit solution for x' > which gives the transition matrix T, is 


C (+) = xO l +) + 


E - H’ + ie V AE 


17T 


W^Qp)^ 


(B13b) 


Comparing equation (B13a) and (B13b) shows that 


U CN X ^ “ "^^PqA-^QP^o^ ( b14 ) 

where A = (AE) *A. Inserting the reciprocal form of equation (B13b) into equation (B14) 
gives 


U CN ~ -W-^pqA^qp 


1 - iff 


1F7“ ^pq^qp 


(B15) 
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APPENDIX C 


PICKUP PROCESSES 

The projection operators previously developed are applied to a concrete example, 
namely, for the case where two kinds of open channels are related through a pickup (or 
stripping) process. Only a single channel of each type is assumed to be open; that is, 
the internal eigenstate in both channels is that in which all composite fragments occupy 
their ground states only. The extension to multiple open channels is straightforward. 
The a fragments consist of the nucleon n and nucleus si. In the (3 set, particle 1 
of nucleus si is bound to n forming the composite particle D, while the remaining 
particles 2, . . ., A form the nucleus si- 1. The following position vectors for part- 
icle i will be employed, where the nucleon n is labeled by i = 0: 

r. position vector of particle i measured from center of mass of total system as 
origin 

from center of mass of nucleus si as origin 
from center of mass of nucleus si - 1 as origin 
771,770 from center of nucleus D as origin 

v=v 1 -v 0 


The spin variables will not be explicitly indicated. The following kinematical relations 
between the position vectors are needed. Here, is the mass of particle i: 


r D 


= (M 0 + M 1 )- 1 


Mji+Mo 




(Cla) 



( 1+m o m a) ? ' 


(Clb) 




M 


A-l 




(Clc) 


The natural set of coordinates in the a channel is [7q, f j, Ig* • • • , T^) anc * in the 0 
channel is (r^, rf, • • • » ?a)' E 1 uations (Cla) to (Clc) give the linear transforma- 
tion between them. The internal state wave functions are 


40 



(C 2a) 


^ao ~ *o(h> ’ ’ ■> ^ a ) 

(p po = ^T>( ri )'l / o(^2’ * * *» O (C2b) 

(The coordinates | . can be treated as independent if a delta function 6^1f.) is included 
in any integration which involves them, and similarly for the £j.) The subspace l 
comprises all functions of the form x(r^4> Q . An orthonormal basis for the subspace is 
the set of all 


x i = X i^o ( C3 > 

where the X. are a complete orthonormal set of functions of Fq: 

E x i( ? o) x f( ? o)“»(V T i>) 

i 

The projector into the space Y^ is given by 


n 


A “ 



(C4) 


Note that is proportional to a unit operator in Fq -space. A complete coordinate 
representation is 

n A - 6 ( F 0 ' p oW*i> • • •> «a) 

Similarly, an orthonormal basis for Y^ is the set of all 

y j - Y ^o < c5 > 

where the completeness of the orthonormal set Yj is expressed by the closure relation 


Z Y j( F D) Y K 7 b) =6 < r D---b) 


The projector Ilg = 1— 


D 


(^D^o in *° * ias tlie re P resen tation 
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n B = < r D - ■ ■ ■■ ?aW(?’) 2- • • • > ?a) < C6 > 

It is convenient to define the function 



where the subscripts indicate the partial set of coordinates that are integrated over in 
forming the scalar product. Since <p Q is a bound state of nucleus si, 

g(Mj) ~ 0 when 1*1 1 » R A (°8a) 

with representing some effective radius for the nucleus si. This implies, in turn, 
that 


g (? D ,v)~0 when J?f| < R D , |r D | - « (C8b) 

If the projection IlgX^ into is substracted from each basis vector x^, the resulting 

set of vectors span the subspace From equations (C3), (C6), and (C7) there 

follows the result 

n B x i = (?2> • • •• ?a) 


where 



Since is a bound state, equation (C8b) shows that 

*i(?d) - 0 when l ? D l-- 

At asymptotic separations in channel a, |Fq| — °°, || j | < R^, and, in view of equa- 
tion (Cla), this in turn implies that |ryJ — also. Thus, the vectors 

x i = (1 - n B> X i 
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which are orthogonal to y'g and span ^pg> have the same asymptotic behavior in 
channel a as the original basis vectors x.. An easy calculation yields 






Hence, the xt are not, in general, orthogonal to one another, and, in fact, need not be 
a linearly independent set. Of course, the subspace ^p B does not depend on the 
choice of basis functions X.(rQ), and the construction of the subspace may be facilitated 
by an expecially suitable choice of these functions. 

The subspace '^gp^ is spanned by the vectors 


y] * d - n A )y, 


and it is easily shown that 


"a?! = h j( 7 oK(«l h) 


where 


h j( r o) =(K 7 0> r i) | Y j( F 0’?1>d(V«i)) 5(1) 

Y i( r 0>«l)' Y j( ? D) 


and 


^d( ? 0’H) 3 *d(*) 


Directly from (C8a) and the bound state nature of follows the property 

h 3 (r 0 ) - ° when |r 0 |-« 

which, in turn, implies that hj vanishes at asymptotic separations in channel /3. Also, 

(yj |yj) - Sy - ^ |hj) 

shows that the yj are not orthogonal or necessarily independent. 
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A somewhat different projection space from f^g was employed by Feshbach (ref. 2) 
in treating the pickup process to illustrate his formalism. This space (hereinafter called 
can k e viewed as the larger space obtained by taking the vector sum of subspaces 
like over the full set of states both bound and continuum, of the system D. 

Then corresponds to all states having the form Z^Fg, ?f^ o , and an orthonormal 

basis for the space can be expressed as the totality of functions 

z mn = Z mn( r D’^K 

where the Z mn constitute a complete set of functions of and rf (or Fq and 
A particular choice would be Z mn = Y m is assumed that the Z mn are 

normalized so that the closure relation has the form 

E Z mn( F 0’«'l) 2 S.n( ? 6>ri) - < F o ' p oM?l ‘ «D 

m, n 

Clearly, the space 'f'-pB contains the entire component of ^ on the vector In addi- 

tion to the open /3 channel, this space includes the open three- particle channel (n, 1, 
si - 1) in which the nucleus sd - 1 is in the state if/ Q and, also, the contribution to the 
channel am arising from the component of cp on i // . Thus, can be used for 

a wider range of conditions than the two-channel case considered here. The projector 


into T^ F g is 

n FB = S | Z mn)( z mn| = 1 r 0 1 f 1 K>)(^o| (C9) 

m, n u 1 1 

The projector that selects the open channels can be taken to be 

P F = n FB + n ApFB (ClOa) 

= n A + n pgpA. (ClOb) 

The space ^ppg is spanned by the vectors 

x!' = N(1 - n FB )x. (Cl la) 

= NX i( r o){^o(^l» * - *’ C\) “ ^( r 0’^l) ^0(^2’ • ’ •’ O] (Cllb) 
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whose scalar products are 


N ' 2 ( x i'K> - 6 ii - | x i( F oM F o-fi)) < cl2 > 

Although the x!’ are not, in general, orthogonal to one another, they will be orthogonal 
if the function g happens to be independent of Fq. Then the scalar product on the right 
of equation (Cl 2) will be proportional to the Kronecker delta by virtue of the orthogonality 
of the functions X^. In the limit that the nucleus sd - 1 is infinitely massive compared 
with D, 


and the x!’ 


*i = Si = r i 

g = g(rj) 

constitute an orthogonal set. A 

• • • ’ ? a) 


m a _i » m d 

function can then be defined by 
- N (?o - s( F lK] 


(Cl 3a) 
(Cl 3b) 


and (Cl la) becomes 

x i'= x i( F oK 

The normalization constant N, which makes ($q|0q) equal to 1, satisfies the relation 

|N|" 2 = 1 - (g|g) 

Under these conditions, the explicit expressions for the projection operators in equa- 
tion (ClOa), are as follows: 


n ApFB = 2 pX*!' 
i 

hx*o 


= i- 


Pp - ^F ^F 

* r 0 r l 


+ 1- 


K)(K 

*$(<-;> I 


M 


[ A-1 » M D 


(Cl 4a) 
(Cl 4b) 


The subspace sd 


FBpA 


is spanned by the vectors 
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z ff 

mn 


= V 1 - n A) z mn 
= M mn[ z mn “ b mn( ? oVo] 


with 


b mn( r o) 


«( 1 ) 


Their scalar products are given by the equation 

^mn^m'n^ ( z mn| z m'n') “ ^mm'^nn' " ( b mn( r o) (hui’n' ( r o)) 


If the basis set Z mn is chosen to be 


Z mn = X m( F oH(«l) 


then 

b mn( ? o) =X m( ? o) c n( F o) 


where 


c n( F o) = (s( ? 0-fl) K( r i)) i(1) 

and the X R are an orthonormal set in the variable | If the c^r^ were independent 
of r 0 (i. e. , if they were constants), it would then follow from the orthogonality of the 
X m (r 0 ) that the b mn ^rp^ were orthogonal for different values of m. Again, in the limit 
of an infinitely massive target for which equations (Cl 3a) and (Cl 3b) hold, the required 
conditions are fulfilled and the c n are constants. In this case, the orthonormal set 
XnO^l) can be chosen so that 

e( F i) * Voft) < cl5 > 

where the constant c Q = [(g|g)]* / ^. Then the orthonormal basis for the subspace 
*FBpA 1S 
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where 


z mn * X m( T o)K ( F l> • • • > F a) M A-1 » M D 


n = - c o*o] 

= \i( F l) *^o "-I, 2, ... 

» 0 in> - » ° = o, i, 2, . . . 


The projector into this space is 


n FBpA = E z mn)( z mn = ‘r. E ^<*i. 

m, n u n 


= 1 ?„ 


+ 1 r,^<^o - V)(V 0 


(016) 


In the last step, use was made of the completeness of the set A^r^J . From equa 
tions (C7), (C9), and (Cl 5), 


| X o^o) ~ c o 1iT Fb|^o) 

can be written. Inserting this form into equation (C16) along with equations (C4) and (C9) 
gives 


n FBpA l N l ( c o n FB ‘ C o 1 ) II A( c o ln FB “ c o 1 ) + 11 FB " c o 2lI FB n A n FB 
After some manipulation, this equation becomes 


n FBpA “ -II A + n FB + l N | 2 ^ - n FB^ n A^ ” n FB^ 


(C17) 


or, alternatively, 

n FBpA = (1 - n A> n FB + I N ! 2 ( c o 1 - "fB^a' 1 - n FB> < C18 > 
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When the form given by equation (C17) is employed in equation (ClOb), it immediately 
gives equation (ClOa), as a glance at equations (C14a) and (Clla) will confirm. The form 
of equation (C18) is the resolution of the operator II into a part operating on *FB 
and a part operating on the subspace orthogonal to it. Equation (Cl 8) illustrates the rela- 
tion given by equation (29). 
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